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Diver  visibility:  why  one  can  not  see  as  far? 


Weilin  Hou,  Alan  D.  Weidemann 

Navel  Research  Laboratory,  Code  7334,  Stcnnis  Space  Center,  MS  39529 

ABSTRACT 

Diver  visibility  has  been  one  of  the  key  research  areas  in  underwater  vision  and  imaging  studies.  Its  applications  also 
extend  into  imaging  system  performance  evaluation  and  prediction,  whieh  is  important  in  MIW  and  ASW  operations. 
These  applications  are  often  associated  with  eoastal  ocean  waters,  and  this  is  generally  translated  directly  into  turbidity 
of  the  water  column.  While  mostly  this  is  the  case,  exceptions  can  lead  to  erroneous  predictions  and  p<  tentially 
significant  consequences.  We  examine  issues  associated  with  sueh  situations,  both  by  model  as  well  as  fiek  data,  in 
order  to  reach  better  estimates  and  to  explore  means  to  compensate  for  sueh  effects,  to  enhance  diver  1  isibility. 
Visibility  data  collected  by  Navy  divers  from  clean  and  relatively  ealm  waters  outside  Pensacola,  during  S  pt  2001 
Gorging  Littoral  Ocean  for  Warfighters  (GLOW)  experiments  suggested  a  eloser  examination  is  warranted,  as  observed 
diver  visibility  measured  at  different  spatial  frequencies  eontradiets  conventional  model  predictions.  Observa  ion  data 
from  two  different  days,  by  different  divers  at  different  depths  were  used.  The  modulation  transfer  of  high  f  equeney 
components  disappears  at  a  level  mueh  higher  than  those  predieted  by  the  human  vision  sensitivity  lev  d  Such 
contradictions  can  be  resolved,  once  the  effeet  of  the  turbulenee  scattering  is  considered  using  a  general  imagin  model. 

Keywords*,  visibility,  scattering,  partieles,  turbulence,  underwater,  MCM,  OTF,  MTF 

1.  BACKGROUND 

Diver  visibility  is  one  of  the  key  underwater  vision  and  imaging  topics  dating  back  to  the  early  days  ol  modem 
oeeanographie  research.  One  of  the  simple  approaches  to  determine  diver  visibility,  or  turbidity  of  the  water,  las  been 
the  use  of  Secchi  disk,  for  the  past  150  yearsfl,  2].  The  idea  was  rather  intuitive,  that  the  disappearing  distan  c  of  the 
dinner  dish  would  give  clue  to  the  turbidity  of  the  water.  Its  deployment  is  simple;  one  lowers  the  traditional  y  white, 
circular  disk  about  30cm  in  diameter,  from  above  the  water  into  the  water  eolumn,  and  determines  the  point  at  which  it 
disappears  from  sight.  The  depth  at  which  the  disk  can  no  longer  be  diseemed  against  the  background  is  defin  d  as  the 
Secehi  disk  depth.  The  physics  behind  the  Seeehi  disk  is  the  absorption  but  mostly  scattering  by  the  wate  and  its 
constituents  to  reduce  the  contrast.  The  classical  work  done  by  Duntley  relates  sueh  quantities  by  means  of  adiativc 
transfer[2],  A  recent  study  revealed  sueh  links  from  image  formation,  or  modulation  transfer  perspeetive[l],  c  pecially 
in  horizontal  situations. 

Image  formation  and  transmission  in  underwater  environments  are  always  hindered  by  the  attenuation  by  watt  *  and  its 
constituents.  Unlike  in  the  atmosphere  where  visibility  ean  be  on  the  order  of  miles,  the  visual  range  in  un  lerwater 
environments  is  rather  limited,  at  maximum  on  the  order  of  tens  of  meters,  even  in  the  clearest  waters.  The  c  jminant 
factor  is  typically  not  due  to  absorption,  unless  visible  channels  near  the  red  arc  used.  Scattering  by  w  iter  and 
particulates  within  the  water  often  poses  the  most  challenges,  causing  photons  stray  from  its  designated  p  th,  thus 
introducing  the  blur  and  loss  of  details  within  the  image.  This  reduction  of  details,  or  high  frequency  components, 
prevents  retrieval  of  information  needed  for  detection  and  more  importantly,  for  clear  identification. 

The  main  cause  of  blurry  images  in  water  is  scattering,  ie  the  photon  changes  path  or  direction.  This  is  describe 
volume  scattering  function  (VSF)  which  is  the  probability  of  photons  being  scattered  into  certain  directions.  F 
waters,  the  VSF  has  a  strong  scattering  peak  in  very  small  forward  angle  In  more  turbid  waters,  the  peak  broa< 
larger  angles,  or  higher  chanee  for  photons  scatter  to  larger  angles,  which  means  more  blur[3,  4].  This  also 
multiple  scattering,  where  the  same  photon  being  scattered  more  than  once  before  reaching  detector.  The 
distribution  is  different  from  that  of  the  single  scattering,  with  even  wider  peak  thus  more  blur. 

Such  effects  can  be  best  described  via  the  point  spread  function  (PSF)  of  the  medium  and  system,  or  simply  put,  the 
system  response  to  a  point  source,  denoted  as  h(x,y).  The  outcome  image  g(x,y)  is  the  result  of  the  simple  convo  ution  of 
h(x,y)  and  original  signal  f(x,y),  adding  the  noise  n(x,y): 
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g(x,y) 


(i) 


=  \\f(x„yi)h(x  -  x^y-y^dxjy,  +  n(x,y) , 

Mathematically,  it  is  easier  to  manipulate  the  above  relationship  in  frequency  domain,  as  convolutions  represert  simple 
multiplications  after  Fourier  transforms, 

G(u,  v)  =  F(u,v)H(u,v)+  N(u,v) .  (2) 


The  Fourier  transformed  PSF  is  the  optical  transfer  function  or  OTF,  H(u,v).  u,  v  are  frequency  componen  t  When 
consider  only  the  magnitude  of  the  transform,  which  is  adequate  in  underwater  incoherent  imaging  applicati  >ns,  it  is 
referred  to  as  modulation  transfer  function  or  MTF.  We  know'  that  the  total  OTF  can  be  modeled  in  simple  forms  with 
multiplicative  components  as  in  all  linear  systems,  while  the  effect  of  water  involving  mostly  particulate  scatt<  ring  can 
be  expressed  in  a  closed  formf  1 , 4,  5] 


D(y/)r 
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and 
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where  c,  b  are  attenuation  and  scattering  coefficients  respectively.  Q{)  is  related  to  mean  scattering  angle.  Ti  e  single 
scattering  albedo  o>o -b/c  is  used  to  describe  the  amount  of  scattering,  while  the  optieal  depth  i=c*r  (r  is  ran  ;e).  The 
angular  spatial  frequency  is  expressed  in  cycles  per  radian. 


Degradation  of  the  image  quality  in  a  scattering  medium  involving  turbulence  has  been  studied  mostly  in  atmosphere. 
These  studies  are  mainly  focused  on  modeling  the  optical  transfer  function,  in  an  effort  to  restore  the  images  obtained, 
such  as  in  air  reconnaissance  or  astronomy  studies  [6,  7].  Little  has  been  done  regarding  the  turbulence  effects  on 
imaging  formation  in  water,  mainly  due  to  the  dominant  particle  scattering  and  strong  attenuation  associated  1  his  is  of 
little  surprise  as  anyone  with  experiences  in  coastal  waters,  especially  those  inside  a  harbor,  or  estuary  a  eas  like 
Mississippi,  experience  first-hand  looks  of  how  visibility  could  quiekly  reduee  to  zero  in  a  matter  of  a  few  1  ;et.  The 
same  applies  to  regions  of  strong  re-suspensions  from  the  bottom,  both  in  coastal  regions  as  well  as  in  the  c  ;ep  sea. 
However,  the  effects  of  turbulence  have  been  postulated  to  have  impacts  over  long  image  transmission  range  [5],  and 
supported  by  light  scattering  measurements  and  simulations[8].  Under  extreme  conditions,  observations  have  be  m  made 
that  involves  targets  with  a  few  feet[9].  The  images  obtained  under  sueh  conditions  are  often  severely  deg  aded  or 
blurred,  on  par  or  more  than  those  caused  by  particles. 


2.  A  DIVER  VISIBILITY  PARADOX 

Diver  visibility  can  be  measured  simply  by  using  Secchi  disk  as  mentioned  above,  or  use  more  sophisticate  targets 
involving  different  spatial  frequencies,  sueh  as  the  USAF-type  resolution  eharts  (Fig.  1).  The  benefits  of  the  later  is  that 
it  provides  convenient  spatial  frequency  modulation  measurements,  which  can  be  used  directly  with  underwater  maging 
models[l]  During  Gorging  Littoral  Ocean  for  Warfighters  (GLOW,  Sept,  2001)  experiments,  measurement  we  *e  made 
to  estimate  diver  visibility  under  different  range  and  optical  conditions.  Some  of  the  data  is  shown  in  Table  1 . 


It  can  be  readily  seen  by  above  visibility  model  (Eqs.  (3)  and  (4))  that  at  the  range  of  disappearance,  the  me  lulation 
transfer  (Hwatcr)  or  relative  contrast  can  be  as  high  as  0.24  for  high  frequency  components  (Fig.2),  which  is  n  direct 
contradiction  of  common  sense  and  those  of  Blackwell’s  eonclusions[10].  In  other  words,  the  imaging  model  prec  lets  that 
the  diver  should  still  be  able  to  see  the  high  frequency  group  (finer  line  pairs),  but  the  diver  reported  otherw  se.  The 
reason  behind  this  contradiction  could  be  the  effects  of  turbulence  as  postulated  above.  We  will  apply  a  ecently 
developed  imaging  model  in  the  following  section  to  solve  the  puzzle. 
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Tabic  1 .  Diver  observed  disappearing  distance  in  feet,  f*  denotes  spatial  frequency  or  line  pairs,  corre  ted  for 
water.  Each  observation  series  were  carried  out  by  different  Navy  divers.  No  vision  issues  are  involved. 


Figure  1  Sample  diver  visibility  panel  at  close  range.  Bottom  right  comer  shows  the  enlarged  version  of  the 
resolution  panel  Picture  was  taken  underwater  during  GLOW  experiment. 


3.  APPLYING  GENERAL  UNDERWATER  IMAGING  EQUATION 

A  general  underwater  imaging  model  has  been  developed,  to  include  the  contributions  of  both  particle  and  tu  bulcncc 
scattering  effects  on  imaging  in  the  underwater  environment[l  1].  For  convenience,  we  will  briefly  outline  the  key 
results  before  apply  the  model  to  the  current  problem. 

From  Fourier  optics,  it  is  commonly  known  that  spatial  coherence  functions  between  optical  fields  of  any  two  pc  ints  can 
be  used  to  describe  the  irradiance  distribution  of  the  source  image  or  object[12,  13]  The  famous  Young  inte  ference 
experiment  can  be  seen  as  a  special  case.  The  modulation  transfer  function  can  be  shown  as  equivalent  to  th  spatial 
correlation  function  on  the  pupil  screen.  For  a  time-varying  correlation  function  under  wide-sense  stationary  conditions. 
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its  ensemble  average  can  be  related  to  the  spatial  phase  structure  function  when  absorption  variation  can  be  neglected, 
such  that 


HA*)  =  r,(3c)  =  exp(-£>,  (x)  /  2) 

(5) 

£>,(*)  =  (0(*)-* Hx  +  Ax))2 

(6) 

where  <f>(x )  is  the  phase  of  the  optical  field  on  pupil  plane.  This  is  under  the  assumption  that  amplitude  structure  docs 
not  change  over  the  small  separation  distance  under  examination.  Other  preconditions  involve  small  perturbations  of 
index  of  refraction,  scale  of  turbulence  compared  large  to  wavelength  and  log-normal  distribution  under  weak 
turbulence  can  be  found  in  detailed  discussions  in  listed  references[13,  14]. 

Following  the  Kolmogorov  model,  for  a  fully  developed  turbulent  flow,  under  the  inertial  regime,  the  powei  spectral 
density  of  index  of  refractions  can  be  expressed  in  the  forms  of  [15] 

(7) 


where  the  superscript  K  denotes  Kolmogorov  spectrum.  K/r)  is  equivalent  of  the  structure  constant  of  the  ndex  of 
refraction  fluctuations,  which  describes  the  intensity  of  the  optical  turbulence  strength  (ie  index  of  r .‘fraction 
fluctuations).  It  is  a  function  of  kinetic  energy  dissipation  rate  and  the  molecular  viscosity  of  the  water  The  abo  e  scalar 
relationship  implies  that  the  turbulence  can  be  considered  statistically  isotropic,  and  homogenous  (wi  le-sensc 
stationary).  Further,  it  is  important  to  remember  that  not  all  turbulent  flows  can  be  described  by  the  above 
spectrum[14],  although  experimental  evidence  suggests  the  above  spectrum  can  be  used  in  underwater  conditions  [15], 

Following  the  approach  by  Fricd[16],  and  relating  the  spectral  power  spectrum  to  the  10R  spatial  distribution!  13],  the 
OTF  of  optical  turbulence  in  underwater  environments  can  be  derived  as[  11] 


where 


OTF  (y/,  r)  =  exp 
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The  key  parameter  that  relates  to  turbulence  structure  is  ro,  or  Fries  parameter[17],  which  is  a  function  of  integn  ed 
turbulent  structure  constant  (Kj  or  Cn2  in  equivalent  atmospheric  terms  ),  range  of  propagation  (r)  and  the  optica 
wavelength  involved.  The  general  underwater  imaging  equation  accounts  for  both  particle  and  turbulence  scattc  ing, 
and  takes  the  form’ 
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It  is  easy  to  see  that  OTF(0,r)=exp(-ar).  It  is  also  worth  mentioning  that  in  practice,  the  OTF  is  often  normalhed  by  its 
DC  component.  This  can  be  important  in  situations  especially  when  automatic  gain  control  is  applied 


Assuming  that  the  turbulence  structures  during  GLOW  is  isotropic  and  homogenous,  and  in  fully  developed  inert  al  stage, 
a  characteristic  seeing  parameter  for  underwater  imaging  can  be  introduced  (Rq),  which  is  equivalent  of  Fried  par  Lmcter 
at  Im  distance  in  waterfl  1],  If  one  can  assume  that  the  discrepancies  between  observation  and  model  is  solely  thi  t  of  the 
effects  by  turbulence,  such  effects  can  be  estimated  by  taking  away  the  effects  of  particle  scattering  from  the  tota 
Blackwell  criterion  (ie  2%).  This  sets  /?r0.004. 


Spatial  frequency  (cyc/rad) 


Figure  2.  Effects  of  turbulence  and  particle  scattering  on  OTF.  Top  to  bottom:  turbulence  scattering, 
particle  scattering,  and  the  combined  (multiplicative)  effects. 


The  effects  of  turbulence  on  underwater  imaging  has  been  discussed  and  demonstrated  previously[9  ,  albeit 
qualitatively.  It  has  long  been  postulated  that  such  effects  will  only  pose  significant  challenge  once  high  n  solution 
imaging  becomes  an  issue,  or  in  extreme  turbulent  environments^].  Using  results  above  (Eq.l  1),  we  can  mo<  el  what 
happened  during  GLOW  experiment  using  sample  values  from  Sept  21,  2001  field  measurements.  The  range  is  set  at 
5m.  The  water  has  a  total  attenuation  coefficient  of  0.3m'1.  If  only  particle  scattering  were  present,  the  modulation 
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transfer  of  high  frequency  components  disappears  at  a  level  following  that  of  Eq.  (4),  and  shown  by  the  middle  (green 
color)  curve  in  Figure  2.  This  value  (0.24)  is  much  higher  than  those  predicted  by  the  human  vision  sensitivity  level [  10]. 
Such  contradictions  can  be  resolved,  once  the  effect  of  the  turbulence  scattering  is  considered  (blue  or  top  curv  e),  and 
properly  included  (red  or  bottom  curve). 

Figure  3  includes  all  observations  shown  in  Table  1 ,  using  particle-scattering  only  data  from  Sept  23  for  clarity.  ( )ne  can 
sec  that  with  only  one  exception  of  Sept  21  observation,  all  relative  contrast  values  are  around  0.02,  in  line  with 
Blackwell  criterion.  Contrasting  to  the  total  contribution  of  both  particle  and  turbulence  scattering,  the  particle-only 
results  clearly  show  high  contrast  values  at  high  spatial  frequencies. 
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Figure  3.  Combined  effects  of  particle  and  turbulence  scattering  in  terms  of  contrast  threshold  from  dive*' 
observations  from  two  different  days  (see  text  for  details). 


3.  CONCLUSION 

This  paper  discusses  the  effects  of  turbulence  on  imaging  outcome  in  natural  environments.  Specific  lly,  we 
demonstrate  that  the  previous  diver  visibility  model  based  only  on  particle  scattering  can  lead  to  erroneous  pn  dictions 
under  certain  conditions.  When  the  accumulative  effects  of  turbulence  scattering  is  included,  by  using  thv  general 
underwater  imaging  equation  which  is  based  on  Kolmogorov  power  spectrum,  we  are  able  to  explain  the  «  bscrved 
discrepancy,  assuming  weak  path  radiance  contribution.  Further  validation  of  the  theory  is  necessary,  especial  y  under 
different  turbulence  conditions  and  particulate  concentration.  The  application  range  of  the  developed  theory  sh<  uld  also 
be  tested,  along  w  ith  different  turbulence  models  and  subregimes. 
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